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A  Cr203-forming  Ni-base  superalloy  and  this  alloy  coated  with  a  Pt-modified  aluminide  coating  were 
exposed  to  Si02  powder  and  cyclically  oxidized  at  950  °C.  The  uncoated  alloy  showed  a  considerable 
amount  of  spallation  and  buckling  whereas  the  Pt-NiAl  coated  alloy  remained  protective  throughout 
hundred  1  h-cycles.  The  interfacial  failure  is  mainly  ascribed  to  the  increased  thermal  strain  by  the  encap¬ 
sulation  of  external  Si02  particulates  in  an  outward-growing  Cr203  layer.  However,  the  particles  were 
not  embedded  in  the  thermally  grown  oxide  of  the  Pt-NiAl  coated  alloy  due  to  the  slow  inward-growing 
characteristics  of  A1203  scales.  The  buckling  of  the  Cr203  scale  with  embedded  Si02  was  analyzed  with 
(1)  a  classical  buckling  criterion  using  the  instantaneous  coefficients  of  thermal  expansion  of  the  con¬ 
stituents,  and  (2)  finite  element  analyses  (FEA)  to  estimate  the  local  interfacial  shear  stresses.  It  turns 
out  that  the  thermal  strain  with  embedded  Si02  is  larger  than  the  experimentally  determined  critical 
thermal  strain  (eb)  explaining  the  buckling  of  the  oxide  scale  observed  in  the  experiment.  The  FEA  results 
demonstrate  that  local  shear  stresses  at  the  metal/oxide  interface  are  significantly  amplified  near  the 
Si02  particles  showing  that  the  buckling  of  oxide  can  be  readily  initiated  especially  in  the  vicinity  of  the 
embedded  particles. 

©  2011  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

When  syngas  derived  from  coal  or  other  practicable  sources  are 
used  in  gas  turbines  or  solid  oxide  fuel  cells  (SOFCs),  inorganic  ash 
can  form  and  accumulate  on  interior  hardware  which  may  cause 
accelerated  corrosion  [1,2].  The  composition  of  the  ash  from  coal 
or  coal-derived  gas  combustion  mostly  contains  oxide  compounds 
such  as  Si02,  A1203,  Fe203,  and  CaO.  The  amounts  of  these  deposit 
particulates  are  expected  to  be  reduced  by  clean-up  processes,  but 
vary  depending  on  the  effectiveness  of  such  processes  [3].  The 
nature  of  the  ash  deposit  depends  upon  the  operating  conditions  of 
the  gas  turbine,  particularly  the  turbine  inlet  temperature,  quality 
of  the  fuel,  and  the  characteristics  of  the  intake  air.  For  example, 
in  tests  performed  using  2-3  ppmw  ash  particulate  in  the  combus¬ 
tion  stream,  more  deposition  and  corrosion  were  observed  at  a  gas 
temperature  of  1260°C  compared  with  1093  °C  [1].  The  observed 
corrosion  characteristics  differed  from  either  of  those  of  type  I 
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and  type  II  hot  corrosion  [4,5],  or  of  calcium-magnesium-alumino- 
silicate  (CMAS)-induced  failure  in  thermal  barrier  coatings  (TBCs) 
[6-10].  The  individual  effects  of  CaO  and  CaS04  on  the  degrada¬ 
tion  of  alloys  and  coatings  in  the  temperature  range  between  750 
and  1150°C  were  reported  elsewhere  [11,12].  In  the  case  of  CaO,  it 
directly  reacts  with  Cr203  and  A1203  scales  forming  Ca-Cr  or  Ca-Al 
ternary  compounds,  which  accelerate  degradation  including  sig¬ 
nificant  spallation  under  thermal  cycling.  Although  similar  adverse 
effects  of  such  deposits  are  expected  to  be  observed  in  syngas- 
fueled  SOFCs  as  well,  especially  for  metallic  interconnects,  it  has 
not  yet  been  reported  in  the  literature. 

Si02  is  one  of  the  major  ash  components  from  coal-  or  syngas- 
fired  combustion.  The  concentration  of  Si02  in  the  ash  varies  from 
nearly  zero  in  the  case  of  some  syngas-fired  gas  turbines  [3]  up  to 
approximately  50wt%  in  the  case  of  pulverized  fuel-based  power 
plants  [13,14]  and  oxyfuel  systems  [14,15].  Even  though  Si02  is  not 
detected  from  the  ash  in  some  applications,  it  still  can  be  introduced 
via  other  sources  such  as  the  use  of  other  fuels  or  ingested  air  con¬ 
taining  Si02  during  operation.  It  is  also  important  to  understand 
the  effect  of  the  possible  ranges  of  impurities  for  safe  operation. 
Although  being  a  major  ash  component,  the  effects  of  external  Si02 
on  the  oxidation  characteristics  of  alloys  have  not  been  studied 
since  the  non-reactive  deposit  constituents  have  been  considered 
as  innocuous  species  in  alloy  degradation  probably  because  of 
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their  high  thermodynamic  stability.  However,  it  is  necessary  to 
determine  whether  the  deposits  are  actually  benign  relative  to  the 
degradation  of  high  temperature  alloys  and  coatings  and,  if  not,  the 
mechanism  of  degradation  induced  by  non-reactive  Si02  needs  to 
be  clearly  addressed.  Through  a  combined  effort  comprises  exper¬ 
imental,  analytical,  and  computational  modeling  approaches,  the 
present  study  explores  detrimental  effects  of  the  external  Si02  par¬ 
ticulates  on  the  interfacial  failure  of  Cr203  oxide  layers. 


2.  Cyclic  oxidation  experiments 

2.1.  Experimental  details 

To  study  the  effect  of  Si02  deposits  on  the  degradation  of  tur¬ 
bine  materials,  laboratory  cyclic  oxidation  tests  with  and  without 
Si02  deposits  were  performed.  The  alloy  used  in  this  study  is  a 
directionally  solidified  (DS)  Ni-base  superalloy,  GTD-1 11,  which  is 
widely  used  as  a  first-  and  second-stage  airfoil  material  in  electric 
power  generation  gas  turbines  [  1 6,1 7].  The  nominal  composition  of 
the  alloy  is  Ni-14Cr-3Al-9.5Co-2.8Ta-3.8W-l.5Mo-4.9Ti  in  wt%. 
In  addition,  the  alloy  coated  with  a  platinum-modified  nickel  alu- 
minide  (Pt-NiAl)  was  studied  to  compare  the  differences  between 
outward-  and  inward-growing  oxide  layers;  GTD-1 1 1  alloy  is  a  pre¬ 
dominantly  outward-growing  Cr203 -former,  whereas  the  Pt-NiAl 
coating  is  a  predominantly  inward-growing  Al203-former.  The 
transport  mechanism  of  cations  and  anions  in  Cr203  and  A1203 
scales  determining  the  growth  directions  are  well  described  in  the 
previous  literature  [18].  Pt-NiAl  coatings  were  prepared  by  Alcoa 
Howmet  Inc.  via  a  low-activity  chemical  vapor  deposition  (CVD) 
diffusion  aluminizing  process  of  Pt-coated  GTD-1 11.  Pt-NiAl  coat¬ 
ings  are  mainly  used  in  the  hottest  sections  of  gas  turbines  and  as 
bond  coats  for  thermal  barrier  coatings  (TBCs).  The  roles  and  failure 
mechanisms  of  these  bond  coats  in  TBC  systems  are  well  described 
in  the  literature  [19-21]. 

The  specimens  were  wet  ground  through  1 200  grit  SiC  papers  to 
establish  controlled  surface  condition  for  all  test  specimens.  Sub¬ 
sequently,  the  specimens  were  grit-blasted  at  207  kPa  (=30psi) 
for  45  min  using  3  p>m  a-Al203  powder.  Pt-NiAl  coated  specimens 
were  subjected  to  the  same  grit-blasting  condition,  but  without 
pre-polishing.  This  provided  surface  roughnesses  (fta)  of  uncoated 
and  Pt-NiAl  coated  specimens  which  were  about  0.3  and  0.5  p,m, 
respectively.  After  grit-blasting,  the  specimens  were  rinsed  with 
water  in  an  ultrasonic  bath.  The  Si02  deposits  (1.0  ±0.1  mg  cm-2) 
were  then  applied  to  the  coupon  specimens  by  spraying  an  aqueous 
solution  containing  the  deposit  onto  all  sides  of  each  specimen.  The 
deposited  specimens  were  vertically  hung  in  a  sample  holder  dur¬ 
ing  oxidation.  The  particle  size  distribution  of  Si02  powder  used  in 
this  study  was  characterized  using  Photo-Correlation  Spectroscopy 
(PCS,  Horiba  LB  550)  as  shown  in  Fig.  1.  The  average  particle  size 
was  0.98  |jim  in  diameter  and  the  particle  size  varied  between  0.2 
and  6.0  p>m.  The  actual  ash  particle  size  range  in  the  power  plant 
falls  between  0.1  and  60  p>m  depending  on  fuel  specifications  and 
filtering  processes.  The  particle  size  and  its  distribution  used  in  the 
present  experiments,  therefore,  will  represent  the  smaller  side  of 
the  practical  range.  The  characterization  of  the  particle  size  and 
its  distribution  is  important  since  it  is  directly  related  to  the  con¬ 
tact  area  of  the  underlying  materials  which  may  affect  the  chemical 
reaction  kinetics  and  it  is  also  related  to  the  possibility  of  encapsu¬ 
lation  in  an  outward-growing  oxide  layer. 

The  cyclic  oxidation  experiments  were  performed  at  950  °C  in 
furnaces  with  controlled  gas  environments.  The  specimens  with 
dimensions  of  15  mm  x  10  mm  x  3  mm  were  cycled  in  and  out 
of  the  furnace  automatically;  45  min  in  the  hot  zone,  15  min  in 
the  cold  zone.  The  exposures  were  interrupted  to  measure  mass 
changes  every  twenty  1  h-cycles.  The  deposits  remaining  unre¬ 


acted  on  the  specimen  were  removed  in  an  ultrasonic  bath  prior 
to  weight  change  measurements.  Then,  the  Si02  was  re-deposited 
for  the  following  20  h  run.  After  hundred  1  h-cycles,  exposures  were 
suspended  and  the  specimens  were  characterized  via  optical  met¬ 
allography  (OM),  X-ray  diffractometry  (XRD)  and  scanning  electron 
microscopy  (SEM)  equipped  with  an  energy-dispersive  X-ray  spec¬ 
troscope  (EDAX). 

2.2.  Specific  mass  changes  and  SEM  observation 

Fig.  2  illustrates  the  area  specific  mass  changes,  Amsp  (mg  cm-2 ), 
as  a  function  of  the  number  of  cycles  for  uncoated  (GO)  and  Pt-NiAl 
coated  (GP)  GTD-1 11  specimens  with  Si02  deposits  exposed  at 
950  °C  for  1 00  cycles  in  dry  air.  The  mass  changes  of  the  same  spec¬ 
imens  without  Si02  deposits  are  also  plotted  for  comparison.  The 
kinetic  curve  clearly  shows  that  the  deposits  have  a  negligible  effect 
on  the  mass  changes  of  Pt-NiAl  coated  GTD-1 1 1  (curves  (1 )  GP  and 
(2)  GP  +  Si02)  whereas  there  is  a  considerable  effect  on  those  of 
uncoated  alloy  (curves  (3)  GO  and  (4)  G0  +  SiO2)  especially  after  40 
cycles.  The  abrupt  drop  of  specific  mass  change  implies  the  possi¬ 
bility  of  spallation.  The  Cr203  on  GO  grows  faster  than  the  A1203 


Number  of  lh-cycles 

Fig.  2.  Specific  mass  changes  versus  number  of  1  h-cycles.  Uncoated  GTD-1 1 1  (GO) 
with  SiC>2  showed  an  indication  of  spallation  whereas  the  effect  of  SiC>2  on  oxidation 
of  Pt-NiAl  coated  GTD-1 11  (GP)  was  negligible. 
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on  GP  because  the  ionic  diffusivity  of  Cr3+  through  the  Cr203  layer 
(outward-growing)  is  much  larger  than  that  of  O2-  across  the  A1203 
layer  (inward-growing). 

To  further  characterize  the  different  kinetic  behaviors  of  the 
oxide  layers,  the  surfaces  and  cross-sections  of  the  specimens  with 
and  without  particles  were  examined  by  SEM.  Fig.  3  shows  sur¬ 
face  SEM  micrographs  for  (a)  GO  specimen  without  Si02  and  (b) 
GO  with  Si02  exposed  at  950  °C  after  100  cycles.  Significant  spal¬ 
lation  of  the  oxide  scale  was  evident  in  the  presence  of  Si02,  as 
shown  in  Fig.  3(b)  at  low  magnification  (50x).  Local  areas  were 
enlarged  to  closely  reveal  (c)  a  spalled  region,  and  (d)  aTGO  retained 
region  at  higher  magnification  (2000 x).  The  micrograph  for  the 
spalled  region,  Fig.  3(c),  was  characterized  via  SEM  and  EDAX  as  a 
polycrystalline  grain  structure  of  the  substrate  alloy  with  fine  Si02 
particles.  In  Fig.  3(d),  it  was  found  that  external  Si02  particles  were 
embedded  in  the  Cr203  in  the  regions  where  the  Cr203  is  retained 
on  the  surface.  This  observation  is  interesting  since  the  encap¬ 
sulation  phenomenon  has  not  been  previously  reported  to  our 
knowledge. 

Cross-section  micrographs  of  the  GO  with  Si02  after  100  cycles 
are  presented  in  Fig.  4  which  shows  significant  amounts  of  scale 
buckling,  Fig.  4(a),  and  crack  opening  at  the  oxide/alloy  interface, 
Fig.  4(b).  Also,  it  is  shown  that  Si02  particles  were  encapsulated 
in  the  Cr203  as  observed  from  the  surface  SEM  micrographs  in 
Fig.  3.  For  comparison,  the  cross-section  micrograph  for  the  same 
specimen  (GO)  exposed  under  the  same  exposure  condition,  but 
without  Si02  deposits,  is  shown  in  Fig.  4(c),  indicating  the  adherent 
oxide/metal  interface  with  a  larger  TGO  thickness. 

When  the  Si02  particles  are  embedded  in  the  Cr203  TGO,  it  is 
expected  that  the  difference  in  the  coefficients  of  thermal  expan¬ 
sion  (CTEs)  between  oxide  and  alloy  ( Aa)  will  greatly  increase  since 
the  mean  CTE  of  crystalline  Si02  (Qf  =  ~0.55  x  10~6I<_1  [22])  is  more 
than  one  order  smaller  than  that  of  Cr203  (of  =  ~7.5  x  10-6  K-1  [22]). 


Estimation  of  the  critical  strain  induced  by  the  Aa  to  initiate  buck¬ 
ling  ( sb )  in  GO  with  embedded  Si02  particles  will  be  discussed  in 
more  detail  in  the  following  section. 

On  the  other  hand,  the  effects  of  the  external  Si02  deposits  on 
the  oxidation  characteristics  of  the  GP  specimens  were  quite  differ¬ 
ent  from  the  previously  described  effects  on  GO.  Fig.  5  shows  surface 
and  cross-section  micrographs  of  GP  specimens  with  and  without 
Si02  deposits  cyclically  exposed  at  950  °C  for  100  cycles.  No  indica¬ 
tion  of  spallation  and  encapsulation  of  the  particles  was  observed, 
as  expected  from  the  kinetic  curves  in  Fig.  2.  The  typical  grain 
boundary  ridges  of  Pt-NiAl  coatings  [23]  were  still  evident  with  a 
small  amount  of  fine  Si02  particles  sitting  on  the  surface  as  shown  in 
Fig.  5(c).  The  observation  was  also  consistent  in  the  cross-sections, 
Fig.  5(d),  showing  that  the  alloy/oxide  interface  is  adherent,  and 
only  a  small  amount  of  fine  Si02  particles  was  observed  at  the  top 
of  the  oxide  scale.  Therefore,  the  oxide  layer  remained  protective, 
similar  to  the  case  of  the  layer  with  no  deposits.  No  major  dif¬ 
ference  was  observed  between  specimens  with  and  without  Si02 
when  comparing  Fig.  5(a)  for  GP  and  Fig.  5(b)  for  GP  with  Si02.  The 
reason  that  the  particles  have  not  been  encapsulated  in  the  TGO  can 
be  explained  by  the  very  slow-  and  inward-growing  characteristics 
of  the  polycrystalline  A1203.  A1203  is  known  to  be  a  primarily  slow 
inward-growing  oxide  because  of  its  dense  structure  with  a  low 
defect  concentration  and  the  smaller  diffusivity  of  Al3+  compared 
with  the  diffusivity  of  O2-  across  the  oxide  layer,  mainly  through 
its  grain  boundaries.  However,  the  A1203  still  grows  outward  to  a 
small  extent  via  the  outward  cationic  transport  through  the  oxide 
lattice.  This  growing  oxide  can  ‘grab’  fine  particles  and/or  lower 
parts  of  larger  particles  during  oxidation.  For  example,  small  Si02 
particles  sitting  on  the  growing  A1203  oxide  surface  are  clearly  seen 
in  Fig.  5(c).  The  difference  in  the  encapsulation  behavior  between 
(a)  a  fast  outward-growing  and  (b)  a  slow  inward-growing  thermal 
oxide  is  illustrated  in  Fig.  6.  It  shows  schematically  that  external 
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Fig.  3.  Surface  SEM  micrographs  of  uncoated  GTD-111  (GO)  with  SiC>2  deposits  exposed  at  950  °C  for  100  cycles  in  air.  Spallation  is  evident  in  (b)  and  (c).  SiC>2  particles  are 
embedded  in  a  thermally  grown  Cr203  layer  as  shown  in  (d). 
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Fig.  4.  Cross-section  SEM  micrographs  of  uncoated  GTD-111  (GO)  cyclically  exposed  at  950  °C  in  air  after  100  cycles  with  and  without  Si02  deposits.  It  is  shown  that  (a) 
buckling  and  (b)  cracking  at  the  metal/oxide  interface  can  occur  in  the  presence  of  Si02,  whereas  (c)  an  adherent  and  protective  Cr203  scale  is  evident  in  the  absence  of  Si02. 


(c)  GP  +  Si02  at  higher  magnificatio 
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Fig.  5.  Surface  and  cross-section  SEM  micrographs  of  Pt-NiAl  coated  GTD-111  specimens  (GP)  with  and  without  Si02  deposits  exposed  at  950  °C  for  100  cycles  in  air  to  show 
that  (b)  neither  significant  spallation  nor  Si02  particle  encapsulation  in  a  TGO  occurred  over  the  whole  surface  and  (c)  grain  boundary  ridges  are  still  evident  indicating  that 
only  a  thin  Al203  layer  is  formed  on  the  surface.  It  is  seen  that  small  Si02  particles  are  still  sitting  on  the  surface.  The  surface  of  GP  without  Si02  is  also  shown  in  (a)  for 
comparison.  A  cross-section  micrograph  of  the  specimen  is  presented  in  (d)  showing  that  there  is  no  indication  of  Si02  encapsulation. 
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(a)  on  a  fast  outward-growing  scale  (b)  on  a  slow  inward-growing  scale 


V 

a  particle  TGO  external 

\ 

external 

particles 

sitting  on  a  particles 

TGO  substrate 

original 

surface 


Fig.  6.  Cartoons  illustrating  that  (a)  external  particles  can  be  encapsulated  in  a  fast  outward-grown  oxide  such  as  Cr2  03  scales,  and  (b)  not  in  a  slow  inward-growing  oxide 
such  as  AI2O3  scales. 


particles  can  be  encapsulated  only  in  a  fast  outward-growing  ther¬ 
mally  grown  oxide. 


3.  Numerical  analyses 

Two  methods  of  analysis  were  adopted  to  understand  the  inter¬ 
facial  failure  of  a  Cr203  scale  in  the  presence  of  embedded  Si02 
particles  with  emphasis  on  buckling.  First,  using  the  classical  buck¬ 
ling  criterion  developed  for  an  ideal  flat  interface,  the  thermal 
strains  arising  from  cooling  were  estimated  and  compared  with  the 
critical  strain  for  buckling  determined  based  on  the  experimental 
observations  of  the  current  study.  Secondly,  Finite  Element  Anal¬ 
ysis  (FEA)  was  utilized  to  investigate  the  local  distributions  of  the 
interfacial  shear  stresses  (rxy)  for  idealized  microstructures  with 
and  without  embedded  Si02  particles.  This  analysis  of  the  failure 
via  embedded  particles  was  developed  based  on  intensive  studies 
by  previous  workers  [24-31  ]. 


3.2.  Determination  of  the  critical  thermal  strain  (sb)  and  the 
critical  temperature  drop  (ATb)for  buckling 


If  it  is  assumed  that  (1)  an  interfacial  crack  opening  occurs  in 
Mode  II,  i.e.,  crack  surfaces  are  displaced  in  sliding  mode  by  shear 
stress  [24],  and  (2)  the  substrate  alloy,  its  TGO,  and  external  partic¬ 
ulates  are  under  elastic  conditions,  one  may  adopt  the  formulation 
for  the  critical  buckling  stress  proposed  by  Timoshenko  and  Gere 

[25]  and  applied  to  cases  of  failure  in  thermal  oxides  by  Evans  et  al. 

[26] .  The  required  stress  to  cause  buckling  is  represented  by: 
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Accordingly,  the  critical  temperature  drop,  A Tb,  for  buckling 
with  an  assumption  of  a  constant  Aa  can  be  estimated  as: 
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Most  previous  work  considered  the  Aa  as  constant  with  respect 
to  temperature  to  calculate  sb  and  A Tb.  In  the  case  of  Al203-forming 
alloys,  for  example,  it  is  reasonable  to  assume  that  the  Aa  is  con¬ 
stant  since  the  CTE  changes,  with  increased  temperature,  of  most 
Al203-forming  alloys  are  almost  parallel  to  that  of  the  A1203  oxide 
layer.  In  other  words,  in  many  cases  of  Al203-forming  alloys,  the 
slopes  of  the  CTE  variations  with  temperature  are  similar  in  alloys 
and  oxides.  As  a  result,  the  CTE  difference  between  the  alloy  and  its 
thermally  grown  A1203  remains  constant  at  different  temperatures. 
However,  in  the  case  of  Cr203 -formers,  the  Aa  between  the  alloy 
and  Cr203  varies  largely  and  non-linearly  in  different  temperature 
ranges.  Therefore,  Aa  should  not  be  assumed  to  be  a  constant,  and 
sb  is  more  appropriate  for  a  buckling  criterion  instead  of  the  crit¬ 
ical  temperature  drop.  This  will  be  further  discussed  in  a  separate 
publication. 

To  analytically  calculate  the  thermal  strains  as  a  function  of  tem¬ 
perature  drop  from  950  °C,  the  existing  CTE  data  collected  from  the 
literature  were  used.  Fig.  7(a)  shows  instantaneous  CTE  data  for 
GTD-111  [32],  Cr203  [22],  and  Si02  [22],  as  well  as  the  calculated 
upper  and  lower  bounds  of  CTE  of  Si02  embedded  Cr203.  Assum¬ 
ing  that  the  Si02/Cr203  composite  is  isotropic  and  the  interfaces 
between  Si02  and  Cr203  are  adherent,  one  can  calculate  the  lower 
(isostress  condition)  and  upper  bound  (isostrain  condition)  of  the 
CTE  for  the  composite  using  the  following  equation  based  on  the 
balance  between  the  total  strain  and  the  sum  of  strains  of  each 
phase: 
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composite 


•  «Si02  +/cr 2O3  ■  <r203 


(4) 


in  which  ob  is  the  critical  buckling  stress;  Eox  is  the  elastic  modulus 
of  the  oxide;  v0x  is  the  Poisson’s  ratio  of  the  oxide;  h  is  the  thickness 
of  the  scale;  and  R  is  the  radius  of  the  interfacial  zone  of  decohesion, 
respectively.  If  the  only  source  of  stress  is  derived  from  differential 
thermal  contraction,  the  critical  thermal  strain  for  buckling  is  given 
by: 

sb  =  p-  =  ATb-  Aa  (2) 

Eox 

where  sb  is  the  critical  strain  to  cause  buckling  and  Aa  is  the  differ¬ 
ence  between  linear  thermal  expansion  coefficients  of  metal  and 
oxide  (=am  -aox).  Insertion  of  Eq.  (1)  into  Eq.  (2)  gives  a  critical 
thermal  strain  to  initiate  buckling  (eb),  which  is: 

•.-AVA.-^)^)’  ,3a, 


in  which  /si02  and  fcr2o3  are  the  volume  fractions  of  Si02  and 
Cr203,  respectively;  n  is  1  for  the  isostrain  condition  and  -1  for 
the  isostress  condition.  Estimation  of  volume  fraction  of  a  specific 
component  from  a  cross-sectional  area  is  typically  accomplished 
by  using  the  volume-area  proportionality  relationship  proposed 
by  Delesse  [33].  It  enables  the  measurement  of  the  volume  fraction 
of  encapsulated  Si02  in  a  TGO  via  a  two-dimensional  pixel  counting 
method  from  a  digitized  map  of  a  representative  micrograph  (see 
Fig.  9(a))  providing  a  value  of  about  0.29.  Based  on  the  collected 
and  calculated  data  given  in  Fig.  7(a),  the  Aa  between  GTD-111 
and  Cr203  with  and  without  embedded  Si02  particles  as  a  func¬ 
tion  of  temperature  drop  from  950  °C  were  obtained  as  shown  in 
Fig.  7(b).  The  area  under  each  curve  in  Fig.  7(b)  corresponds  to  a 
thermal  strain,  which  can  be  integrated  as  a  middle  Riemann  sum. 
Fig.  8  plots  these  calculated  thermal  strains  as  a  function  of  tem¬ 
perature  drop  from  the  oxidation  temperature  of  950  °C  with  and 
without  the  embedded  Si02  particulates  at  every  50 1<  interval. 
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Temperature  drop  from  950°C,  AT  /  K 

Fig.  7.  (a)  Instantaneous  CTEs  of  GTD-111,  Cr203,  and  Si02  as  well  as  calculated 
upper-  and  lower-bounds  of  CTEs  of  a  Si02-embedded  Cr203  composite  and  (b)  CTE 
differences  with  respect  to  temperature  drop  from  950  °C  between  GTD-111  and 
oxide  scales  with  and  without  embedded  SiC>2.  Note  that  Aa  with  embedded  Si02 
particulates  increased  because  of  the  low  CTE  of  SiC>2. 
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Fig.  8.  Calculated  thermal  strains  as  a  function  of  temperature  drop  from  950  °C 
using  data  in  Fig.  7  along  with  the  critical  thermal  strain  (sb  =  0.80%  elongation) 
obtained  from  Eq.  (3a)  using  measured  values  of  parameters  from  a  cross-section 
micrograph. 


Table  1 

Elastic  properties  used  in  the  present  study. 


Material 

Elastic  modulus  (at  600 °C,  MPa) 

Poisson’s  ratio 

GTD-1 1 1 

140.0  [27] 

0.29  [27] 

Cr2  03 

260.0  [24] 

0.30  [24] 

Si02 

65.7  [24] 

0.20  [24] 

Additionally,  the  critical  thermal  strain  to  initiate  buckling  in  the 
present  alloy/oxide  system  can  be  determined  using  Eq.  (3a)  with 
measured  values  of  the  thickness  of  the  scale  (ft  =  4.1  ptm  on  aver¬ 
age)  and  a  radius  of  a  buckled  zone  ( ft  =  53.0  ptm)  from  the  observed 
microstructure  such  as  Fig.  4(a).  When  the  Poisson’s  ratio  of  the 
oxide  scale  is  0.29,  the  equation  gives  a  value  of  sb  as  0.80%.  This 
means  that  the  oxide  scale  will  buckle  when  the  thermal  strain 
caused  by  the  temperature  drop  is  larger  than  0.80%.  The  sb  calcu¬ 
lated  from  a  buckled  microstructure  is  designated  as  a  horizontal 
dashed  line  in  Fig.  8.  It  shows  that  buckling  is  expected  to  occur, 
with  embedded  Si02,  in  the  temperature  drop  range  476-740 K 
(indicated  by  the  two  vertical  dotted  lines)  as  observed  in  Fig.  4(a). 
The  temperature  drop  to  initiate  buckling  is  referred  as  the  critical 
temperature  drop  for  buckling,  A  Tb.  However,  in  the  absence  of  the 
Si02,  the  calculation  predicts  that  the  scale  would  not  buckle  up  to 
A T=  900  K  consistent  with  the  observation  in  Fig.  4(c),  which  shows 
the  adherent  interface  between  oxide  scale  and  substrate  alloy. 
Note  that  the  value  of  sb  could  vary  depending  on  the  diameter  of 
the  observed  buckled  region,  ft ,  and  the  local  variation  of  the  scale 
thickness.  Sampling  of  the  several  cross-sections  indicated  that  the 
diameter  of  the  buckled  zone  varies  between  44.7  and  76.1  pim; 
thus,  it  is  not  unreasonable  to  assume  the  measured  value  from  the 
cross-section  (ft  =  53.0  pim  and  ft  =  4.1  pun  on  average  over  a  long 
range)  can  be  taken  as  a  representative  of  the  buckled  zone.  The 
elastic  properties  used  in  the  calculation  are  summarized  in  Table  1 . 
Although  the  elastic  properties  of  directionally  solidified  GTD-1 1 1 
will  vary  along  different  loading  directions  [34],  the  averaged  value 
over  all  loading  directions  was  assumed  as  a  representative  of  the 
alloy. 

Based  on  the  calculational  analyses,  it  can  be  additionally  pre¬ 
dicted  that  the  dramatic  effect  of  the  encapsulated  particles  can 
be  reduced  in  the  case  of  the  external  particles  with  higher  CTEs 
than  those  of  Si02  (a  =  ~0.55  x  10_6I<_1),  such  as  Fe203  deposits 
(of  =  ~13.8  x  10_6I<_1  [22]),  which  are  closer  to  those  of  the  sub¬ 
strate  alloy  since  it  will  decrease  the  CTE  difference  between  metal 
and  oxide  leading  to  smaller  thermal  strains. 

Although  this  classical  buckling  criterion  can  be  used  to  evaluate 
this  type  of  interfacial  failure  caused  by  the  encapsulated  particles 
in  a  thermal  oxide,  several  issues  must  be  resolved.  For  example, 
important  aspects  that  need  to  be  considered  are  the  amount  of 
particulates  encapsulated  in  a  thermal  oxide  and  the  cause  of  this 
encapsulation.  It  seems  that  the  volume  fraction  of  embedded  par¬ 
ticulates  varies  as  a  function  of  growth  rate  of  a  thermal  oxide. 
If  this  is  the  case,  oxidation  temperature,  time,  and  kinetic  con¬ 
stants  associated  with  defect  structure  of  the  oxide  will  be  related 
to  the  volume  fraction  of  embedded  particles.  The  ratio  of  par¬ 
ticulate  size  to  oxide  thickness  (r/ft)  will  be  important  as  well. 
Studies  to  determine  the  factors  that  control  the  volume  fraction 
of  embedded  particles  in  a  thermal  oxide  during  growth  may  be 
interesting. 

3.2.  Calculation  of  the  local  shear  stresses  at  metal/oxide 
interface  via  FEA 

To  further  understand  the  localized  effect  of  embedded  Si02 
particles  in  a  Cr203  scale  on  interfacial  failure  of  the  oxide  scale, 
the  changes  in  stress  distributions  during  cooling  were  calculated 
using  two-dimensional  finite  element  analysis  (FEA).  The  calcu- 
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lated  results  via  the  FEA  will  be  discussed  in  terms  of  the  shear 
stresses  (rxy)  at  the  metal/oxide  interface  that  should  play  an 
important  role  on  decohesion  of  the  oxide  scale  from  the  sub¬ 
strate  required  for  buckling  failure.  The  discussion  starts  with  a 
description  of  a  method  to  create  idealized  microstructures  from 
experimental  observation  to  use  as  an  initial  microstructure  for 
the  simulation.  Subsequently,  details  for  the  cooling  simulation  are 
described.  The  interfacial  shear  stresses  calculated  from  the  simu¬ 
lation  will  be  discussed  to  elucidate  the  role  of  embedded  particles 
on  decohesion  of  the  oxide  scale. 

3.2.1.  Generation  of  an  idealized  microstructure  from  a  real 
microstructure 

To  perform  a  more  realistic  calculation,  one  of  the  cross- 
sectional  microstructures  representing  interfacial  delamination 
in  Fig.  4(b)  was  selected  and  digitized  for  the  calculation.  All 
the  major  interfaces  among  embedded  Si02  particulates,  a  Cr203 
scale,  and  a  GTD-111  substrate  were  traced  to  create  a  simplified 
microstructure  as  in  Fig.  9(a).  For  the  generation  of  the  simplified 
microstructure,  it  was  assumed  that:  (1)  there  is  no  pre-existing 
crack  at  the  oxide/metal  interface  and,  therefore,  the  interface 
is  adherent  in  the  digitized  microstructure  and  (2)  only  major 
constituents  (substrate  alloy,  Cr203  TGO,  and  large  embedded 
Si02  particulates)  were  included  in  the  simplified  microstruc¬ 
ture,  i.e.,  all  other  minor  constituents  in  the  oxide  and  substrate 
were  neglected.  The  traced  microstructure  was  implanted  into  an 
FEM  mesh  via  the  OOF  package  (the  Object-Oriented  Finite  ele¬ 
ment  code  developed  by  the  National  Institute  of  Standards  and 
Technology,  Gaithersburg,  MD  [35])  using  an  adaptive  mesh  gen¬ 
eration  algorithm.  The  resulting  mesh  structures  created  by  the 


described  protocol  are  shown  in  Fig.  9(a)  and  (b)  representing 
the  idealized  microstructures  with  and  without  embedded  Si02, 
respectively.  The  microstructure  without  embedded  particles  in 
Fig.  9(b)  was  generated  by  eliminating  the  embedded  particles  from 
the  microstructure  in  Fig.  9(a). 

3.2.2.  Local  shear  stress  distribution  along  metal/oxide  interface 
As  previously  noted,  the  residual  shear  stress  distribution  (rxy) 
arising  from  temperature  changes  is  an  important  measure  to 
describe  ‘decohesion’  of  an  oxide  layer  from  the  substrate  alloy. 
The  rxy  distributions  were  calculated  through  a  cooling  simulation 
via  a  commercial  FEA  package,  ABAQUS  (Abaqus  Inc.,  Providence, 
RI)  [36].  During  the  calculation,  it  was  assumed  that  (1)  all  the 
phases  are  homogeneous,  isotropic,  and  ideally  elastic;  and  (2) 
stress-relieving  mechanisms,  such  as  creep,  phase  transformation, 
crack  formation,  are  not  operating  during  the  simulation.  It  was  also 
assumed  that  (3)  the  thermo-mechanical  parameters  for  the  con¬ 
stituent  phases  do  not  vary  with  respect  to  temperature  except  for 
the  CTEs  which  instantaneously  vary  with  temperature.  The  param¬ 
eters  used  in  the  simulation  were  the  same  as  previously  used  in 
Section  3.1.  A  mixture  of  fixed  and  free  boundary  conditions  (BC) 
was  used  for  the  displacement  of  each  element  along  x  (horizon¬ 
tal)  and  y  (vertical)  directions;  the  free  BCs  were  assigned  to  the 
all  boundary  elements  but  the  displacements  of  elements  along 
±x  on  the  left  and  ±y  on  the  bottom  sides  were  set  to  be  fixed. 
To  exclude  the  computational  artifact  for  the  free  boundary  con¬ 
dition  arising  from  a  relatively  small  computational  section,  the 
meshed  microstructure  in  Fig.  9(a)  was  extended  (ten  times  along 
the  +x  direction  using  mirror  images,  and  the  four  times  along  the 
-y  direction  by  enlarging  the  substrate  region)  as  given  in  Fig.  9(c). 


(a)  with  Si02 


alloy/oxide  interface  is 
adherent 


substrate  (GTD-111) 


(b)  without  Si02 

Cr203 

alloy/oxide  interface  is 
adherent 

substrate  (GTD-111) 

,  5mm  j 


(c)  10x4  fold  extended  microstructure  in  a  2D  FEM  mesh. 


Fig.  9.  (a)  An  idealized  microstructure  using  a  real  microstructure  in  Fig.  4(b),  (b)  a  hypothetical  idealized  microstructure  but  with  no  embedded  Si02,  and  (c)  a  10  x  4  fold 
extended  microstructure  implanted  in  a  2D  FEM  mesh. 
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The  first  section  (out  of  ten  repetitive  sections)  of  the  extended 
domain  at  the  top-left  corner  was  then  selected  as  representative 
of  the  simulation  results  as  indicated  with  a  rectangle  in  Fig.  9(c). 
Preliminary  calculations  (not  shown  in  the  paper)  confirmed  that 
this  pseudo-periodic  boundary  condition  effectively  eliminates  the 
spurious  computations.  The  plane  strain  condition  was  used  assum¬ 
ing  that  the  normal  strain  szz  and  the  shear  strains  7 xz  and  7yz 
are  negligibly  small  compared  with  the  cross-sectional  strains.  For 
comparison,  the  same  extension  and  subsequent  calculation  was 
repeated  but  with  a  hypothetical  microstructure  without  embed¬ 
ded  Si02  in  Fig.  9(b).  The  cooling  simulation  extends  from  950  °C 
down  to  20  °C  and  the  stress  distributions  were  monitored  at  every 
50  K  of  the  temperature  drop. 

Fig.  10  shows  the  absolute  values  of  rxy ,  \rxy\,  along  the 
metal/oxide  interface  for  the  microstructures  with  (in  red)  and 
without  (in  blue)  Si02  particles  induced  by  a  temperature  drop  of 
AT  =  -930  K,  i.e.,  at  T=  20  °C.  Since  the  sign  of  rxy  simply  represents 
the  direction  of  the  shear  stress,  their  absolute  values  were  plot¬ 
ted.  The  plot  shows  that  |rxy|  scatters  in  the  ranges  of  ±2.5  and 
±1.2  GPa  for  the  microstructures  with  and  without  Si02  particles, 
respectively.  The  average  values  of  |rxy|  are  475.3  (with  Si02)  and 
322.2  (without  Si02)MPa.  The  result  implies  that  the  presence  of 
embedded  Si02  particles  in  a  Cr203  scale  will  increase  |rxy|  and, 
consequently  lead  to  a  higher  probability  of  the  oxide  decohesion 
from  substrate.  Such  detrimental  impact  by  the  encapsulated  par¬ 
ticles  can  also  be  clarified  by  the  position  of  |rxy|.  The  equilibrated 
microstructure  after  cooling  is  also  shown  on  top  of  the  plot  to 
qualitatively  show  the  dependency  of  the  location  of  the  embed- 


y 


(b)  without  Si02 


Fig.  10.  Calculated  local  thermal  shear  stresses  ( rxy )  along  the  metal/oxide  inter¬ 
faces  after  cooling  (AT=-930K  from  950 °C)  obtained  from  the  equilibrated 
microstructures  (a)  with  and  (b)  without  embedded  SiC>2  particles.  It  is  shown  that 
the  local  shear  stresses  are  significantly  amplified  in  the  presence  of  the  embedded 
Si02,  especially  when  the  particles  are  located  in  the  vicinity  of  the  interface. 


ded  particles  on  the  interfacial  shear  stresses.  It  can  be  seen  that 
|rxy|  at  the  interface  tends  to  be  greatly  amplified  by  about  50% 
on  average  and  locally  up  to  about  300%  when  the  Si02  particu¬ 
lates  are  located  near  the  interface.  Also,  the  resulting  shear  stresses 
seem  to  become  even  larger  near  larger  particles  as  indicated  with 
arrows  in  the  figure.  It  infers  that,  in  the  presence  of  embedded 
Si02  particles,  decohesion  of  the  oxide  scale  can  occur  even  when 
the  average  shear  stress  over  the  entire  interface  is  smaller  than 
the  critical  shear  stress  for  decohesion  because  the  amplified  local 
shear  stresses  near  the  particles  can  be  greater  than  the  critical 
shear  stress.  This  amplification  of  the  shear  stress  will  be  impor¬ 
tant  to  determine  whether  the  interfacial  failure  occurs  in  mode  I 
or  II  when  the  normal  thermal  strain  (£xx)  during  cooling  is  large 
enough  for  both  mechanisms.  If  the  amplification  by  the  presence  of 
Si02  is  critical  to  cause  decohesion,  the  failure  mode  II,  i.e.,  buckling 
of  the  oxide  scale  after  the  preceding  decohesion  at  the  metal/oxide 
interface,  will  operate  preferentially.  Otherwise,  the  failure  mode 
I,  i.e.,  wedging  crack  through  the  oxide  scale  and  the  subsequent 
decohesion,  will  be  predominant. 

Although  the  critical  interfacial  shear  strength  between  Cr203 
and  the  GTD  111  alloy  is  not  currently  available,  it  is  expected  that 
the  amplification  of  the  shear  stress  can  be  critical  to  initiate  crack 
formation  at  the  interface.  For  example,  for  Cr203-forming  Crofer 
22  APU  which  is  an  Fe-22Cr-0.5Mn  alloy  widely  used  as  a  metal¬ 
lic  interconnect  for  SOFC  applications,  the  critical  interfacial  shear 
strength  was  estimated  as  442.5  (STD  =  1 5.2 )  MPa  by  Sun  et  al.  when 
the  Cr203  thickness  is  in  the  range  2.06-2.98  p,m  [37].  Assuming 
that  the  critical  shear  strength  of  Cr203/Crofer  is  not  much  differ¬ 
ent  from  that  of  Cr203/GTD-1 1 1  with  the  TGO  thickness  of  4.1  p,m, 
the  calculated  local  shear  stresses  can  be  roughly  evaluated  to  see 
whether  the  amplified  local  shear  stress  would  be  sufficient  to 
initiate  decohesion.  Recalling  that  the  average  values  of  |rxy|  are 
475.3  (with  Si02)  and  322.2  (without  Si02)  MPa  from  the  FEA,  it  is 
inferred  that  the  presence  of  Si02  can  be  critical  to  cause  interfacial 
decohesion.  For  a  more  accurate  evaluation,  the  shear  stress  distri¬ 
bution  tangential  to  the  local  interface  and  the  shear  strength  for 
this  particular  interface  will  be  required.  It  should  be  noted  that  the 
maximum  shear  stress  (rmax)  at  the  edge  on  the  right  side  of  the 
interface  decreases  as  the  lateral  dimension  of  the  finite  specimen 
increases  as  predicted  by  Tien  and  Davidson;  rmaxocl/L,  where  L 
is  the  width  of  the  specimen  [31].  Also,  the  shear  stress  along  the 
interface  decreases  from  the  edge  to  the  center  of  the  specimen.  The 
|rxy  |  distributions  near  the  center  region,  i.e.,  at  the  squared  region 
in  Fig.  9(c),  represent  minimal  effects  of  particle  encapsulation  in 
the  entire  specimen.  These  Si02  particle  effects  will  be  even  more 
amplified  near  the  edge  region  of  the  specimen  especially  when  the 
specimen  size  is  small. 


4.  Conclusion 

The  present  paper  describes  a  new  interfacial  failure  mode 
induced  by  the  presence  of  chemically  inert  external  Si02  parti¬ 
cles  that  can  be  accumulated  on  interior  hardware  for  syngas-fired 
power  generation.  Accelerated  cyclic  oxidation  tests  were  per¬ 
formed  to  provide  examples  of  the  failure.  The  experimental  results 
showed  that  the  Si02  caused  ( 1 )  spallation  and  buckling  on  a  Cr2 05- 
former  (GTD-1 1 1 )  although  the  amount  of  the  attack  is  much  less 
significant  than  that  induced  by  CaO,  and  (2)  a  negligible  effect  on 
an  Al203-former  (Pt-NiAl).  The  difference  was  explained  by  the 
rapidly  outward-growing  and  the  slowly  inward-growing  char¬ 
acteristics  of  Cr203  and  A1203  layers,  respectively.  The  buckling 
failure  of  the  Si02  embedded  Cr203  TGO  on  GTD-1 11  was  analyti¬ 
cally  accessed  by  a  modified  version  of  a  classical  buckling  criterion. 
The  analyses  showed  that  the  CTE  differences  (Aa)  and  the  result¬ 
ing  thermal  strains  (£th)  between  oxide  and  substrate  alloy  can 
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be  greatly  increased  via  encapsulation  of  external  Si02.  The  cal¬ 
culated  thermal  strains  were  compared  with  the  critical  strain  for 
buckling  ( sb )  obtained  from  the  observed  microstructure  leading 
to  a  prediction  of  the  critical  temperature  drop  for  buckling  (A  Tb) 
as  476-740 1<  in  the  presence  of  embedded  Si02.  This  detrimental 
effect  of  Si02  particles  was  further  investigated  by  FEA  simulation. 
It  is  concluded  from  the  simulation  that  the  localized  shear  stresses 
(r Xy)  at  the  metal/oxide  interface  will  be  significantly  increased 
especially  in  the  vicinity  of  encapsulated  particles  that  can  lead 
to  accelerated  decohesion  which  is  a  pre-requisite  condition  for 
buckling. 

In  summary,  based  on  the  experimental  observations  and  cal- 
culational  analyses,  the  interfacial  failure  in  the  presence  of  Si02 
predominantly  occurs  as  follows: 

(1)  The  external  Si02  particles  are  encapsulated  in  the  outward¬ 
growing  Cr203  scale  during  oxidation. 

(2)  The  overall  CTEs  of  the  oxide  composite  are  decreased  due  to  the 
embedded  particles  with  low  CTEs  such  that  the  CTE  differences 
between  metal  and  oxide,  and  therefore  the  resulting  thermal 
strains  are  increased. 

(3)  The  increased  thermal  strains  lead  to  amplification  of  the  shear 
stresses  (rxy)  at  the  metal/oxide  interface  that  may  result  in 
accelerated  decohesion. 

(4)  When  the  increased  thermal  strains  exceed  the  critical  strain 
for  buckling  ( sb ),  scale  buckling  occurs. 

Although  the  alloy  systems  and  temperature  used  in  the  present 
study  were  taken  from  syngas-fueled  gas  turbines  to  demonstrate 
the  effect  of  chemically  inert  external  oxide  particles,  the  same 
mechanism  can  be  applied  to  other  particle-associated  high  tem¬ 
perature  applications,  such  as  syngas-fueled  SOFCs  or  coal-fired 
oxyfuel  systems. 
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